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Abstract 

The behavior of electrolyte solutions close to a charged surface is studied theoreti- 
cally. A modified Poisson-Boltzmann equation which takes into account the volume 
excluded by the ions in addition to the electrostatic interactions is presented. In 
a formal lattice gas formalism the modified Poisson-Boltzmann equation can be 
obtained from a mean-field approximation of the partition function. In an alter- 
native phenomenological approach, the same equation can be derived by including 
the entropy of the solvent molecules in the free energy. In order to visualize the 
effect of the steric repulsion, a simple case of a single, highly charged, flat surface is 
discussed. This situation resembles recent adsorption experiments of large ions onto 
a charged monolayer. A simple criterion for the importance of the steric effects is 
expressed in terms of the surface charge density and the size of the ions. It is shown 
that when these effects are important a saturated layer is formed near the surface. 
A modified Grahame equation relating the ion concentration at the surface to the 
surface charge density is obtained. 

Key words: Electrolyte solutions; Ion adsorption; Fluid/Fluid interfaces; Stern 
Layer. 
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1 Introduction 



The Poisson-Boltzmann equation is the main tool for studying the behavior 
of ionic solutions [1-4]. Its main advantages are its simplicity, which allows 
for analytical solutions in simple cases, and its surprisingly good agreement 
with experiments. Aqueous solutions of small and macro-ions are of interest 
from an industrial point of view {e.g., colloidal suspensions [5]) and as an 
important component in biological systems {e.g., DNA, charged membranes). 
Therefore, the Poisson-Boltzmann approach was applied to many situations. 
Of particular interest are: (i) adsorption of ions to fiat surfaces [6,7]; (ii) ion 
distribution around a charge cylinder and the so-called Manning condensation 
[8-12]; and (iii) ion distribution around a charged sphere and the so-called 
charge renormalization [10,13]. 

In its simpler form, the linearized Poisson-Boltzmann equation leads to the 
Debye-Hiickel expression, thus providing a simple description of screening ef- 
fects in terms of the Debye-Hiickel screening length [14]. The success of the 
Poisson-Boltzmann approach is quite impressive in view of the various ap- 
proximations which are included in its derivation: it is a mean field approach 
which totally neglects correlations and all specific (non-electrostatic) interac- 
tions between the ions including the ionic finite size. 

Nevertheless, despite its success in describing a wide range of systems, it has 
been known for a long time to have some limitations in several cases: (i) the 
phase transition of electrolyte solutions [15], (ii) the adsorption of charged ions 
to highly charged surfaces [16-19], and (iii) the attractive interactions which 
can be observed between equally charged surfaces in the presence of asymmet- 
ric electrolytes [20,21]. Consequently, there have been numerous attempts to 
improve upon the standard Poisson-Boltzmann approach [22-28]. 

In this study we focus on the second case, where highly charged surfaces attract 
a large amount of free ions from the solution. At high ion densities achieved 
close to the surface, short range ion-ion interactions become comparable to 
the Coulomb interaction and they can no longer be neglected. In particular, 
the ion density is bounded by the maximum value which is obtained when 
the ions are closely packed. Recently, Cuvillier et al. [18] have provided an 
experimental setup which clearly demonstrates this situation. In their exper- 
iments large tungstic acid ions (diameter of about lOA) were adsorbed onto 
a charged monolayer spread at the air/water interface (Fig. 1). The insoluble 
(Langmuir) monolayer consists of charged amphiphilic molecules having two 
moieties: a charged head group favoring the water side of the interface, and 
a hydrocarbon tail favoring the air side. The surface charge density can be 
controlled continuously by changing the monolayer density through a lateral 
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surface pressure. In the experiments, a large discrepancy was found between 
the measured ion concentration near the surface and the high values antici- 
pated by the Poisson-Boltzmann approach. 

Our aim in this study is to include the finite size of the ions in the Poisson- 
Boltzmann approach and study how the ion distribution close to charged sur- 
faces is affected. The standard way of including the finite size of the ions 
in the Poisson-Boltzmann approach is to define a narrow layer close to the 
surface as impenetrable to the ions. This layer is usually referred to as the 
Stern layer [29] and its width is equal to the ion radius. Outside this layer 
the regular Poisson-Boltzmann equation is implemented. In our approach, a 
modified Poisson-Boltzmann equation is derived where the steric forces lead 
to saturation of the ion density at high potentials. This way, the width of the 
saturated layer depends also on the surface charge and is not limited a-priori 
to only one counterion layer. 

In the next section the modified PB equation is derived for different combina- 
tions of ion valencies. Two derivations are presented: first, a systematic path 
integral approach where the ions are put on a discrete lattice, and second, a 
phenomenological free energy formulation where the excluded volume effect is 
added through the entropy of the solvent. This equation is then implemented 
in Sec. 3 to study the adsorption of large ions to flat surfaces. 



2 The Modified Poisson-Boltzmann Equation 

2.1 Lattice Gas Formulation 

Consider an aqueous solution of charged ions. For simplicity we will assume 
that both co-ions and counter-ions have the same size a. This assumption 
can be justifled when all the surface charges are of the same sign, since only 
counterions are then attracted to the surface and reach high charge densities. 
Another simplification is that we do not distinguish between the counterions 
which dissociate from the charged surfaces and the ones originating from the 
added salt. 

For the valencies of the ions we will consider three cases: (i) a symmetric z:z 
electrolyte, (ii) an asymmetric l:z electrolyte, and (iii) 2;-valent counterions 
without additional salt. The different cases will be used to study the applica- 
tion of the modified equation in different physical systems. 

For a symmetric z:z electrolyte the solution contains two charge carriers, with 
charges equal to ±2;e. In order to derive the free energy we will use a discrete 
lattice gas formulation. In this approach, the charge carriers are placed on 
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a three dimensional cubic lattice where the dimensions of a single cell are 
ax ax a (Fig. 2). Thus, by dividing space into discrete cells (lattice sites) and 
limiting the occupation of each cell to a single ion we introduce a short range 
repulsion between the ions. The size of a cell represents the volume of an ion 
up to a numerical pref actor. 

In order to describe the occupation of cells by ions we assign to each cell j, 
which is located at rj, a spin-like variable Sj. This variable can have one of 
three values: Sj = if the cell is empty (occupied by a water molecule), and 
Sj = ±1 according to the sign of the ion that occupies the cell. The partition 
function of the system can now be written in the form [30] 



Sj=0,±l \ ^ j J 

The first term in the exponent is the electrostatic energy, where fc(r) = l/£|r| 
is the Coulomb interaction. The second term is the chemical potential term 
where /ij — ii+ for the positive ions {sj — +1) and fj,j — /i-F for the negative 
ions {sj = — 1). 

The charge density operator can be expressed in terms of the spin variables 
as 



Pc(r) = E^eSj-5(''-ri) (2) 
j 

It is then useful to introduce the density field Pc(i') and the conjugate field 
<^c(r) through the identity 



= J Vpc VlPc exp^i/3 j dr pc^'c - if^ze^Sj^pcirj)^ (3) 

where / Vpc is a functional integral over the values of pc at all space points 
r. It can viewed as the continuum limit of multiple integrals over the values 
of Pc at different points in space (see, e.g., ref. [31]): 

n / dpe(r,) ^ / Vpc (4) 



Using the above identity is equivalent to the Hubbard-Stratonovitch transfor- 
mation (see, e.g., ref. [32]) and leads to 
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j Vpc Vifc exp^-^ j dr dr'pc(r)t;c(r - r')pc(r') 
+ i/3 y dr pc{v)ipc{v)^ 
X exp(-i/5^e^Sj-93c(rj) + ^/5/ijs|) (5) 

Recall that are the discrete coordinates of the lattice sites while r,r' are 
continuous spatial coordinates. It is now possible to trace over the allowed 
values of the spin-like variables, Sj = 0, ±1 (last line in eq. 5). In the continuum 
limit, where physical properties vary on length scales much larger than the size 
of a single site, the sum over the lattice sites can be replaced by a continuous 
integral over space and the partition function simplifies to 

Z^jvpc Vcpc exp(^-^ I dr drVc(r)^;c(r - r')pc(r') 
+ iP J dv pe(r)<^c(r) 

+ J_ y dr ln{l + e/^'^+-^^^^v'c(r) ^ ^i3n-+izpe^c{r)y^ (6) 

Since the exponential is quadratic in pc, its functional integral can be per- 
formed. 

z = y exp(^ - ^ y dr iv<^er 

_l_ J_ y dr ln{l -I- e'^'^+~*^^^'^'=*^''^ -l-e^'*""'"*^'^'^'^'^^'")}^ (7) 



The chemical potentials p± are related to the total number of positive and 
negative ions in the solutions through 



1 dZ I 1 r Q|3^i±l'^z|3eipc{r) 

^ ~ Z d{l3p±) J 1 + e^M+-jz/3e</5c(r) _|_ e^/i-+iz^e<^c(r) 

Where (O) denotes the grand canonical average of the operator O. 

In the bulk, the total number of positive and negative ions is equal, = 
N_ — N/2. It is useful to define the volume fraction occupied by both the co- 

and counter-ions as 0o = Na^/V = 2cba^ where V is the total volume and 
Cf, is the bulk concentration of the electrolyte. In the thermodynamic limit 
N,V ^ oo while C(, and 0o remain finite. Using eq. 8 the chemical potentials 
can be expressed in terms of 
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2 1 - 00 ^ ^ 



In the mean field approximation, the partition function is approximated by 
the value of the functional integral at its saddle point ijj{r) = iip^. The free 
energy of the system is then given by 



= [ dr - ir [dr ln|l + --^cosh[/3;2eV^(r)]| (10) 

OTT J J [1 — 00 J 

where '0(r) satisfies the modified Poisson-Boltzmann equation for a symmetric 
z:z electrolyte [19]: 



2 , STTze Cbsinh(;2/5e'0) 



£ 1 — 00 + 00 cosh(2;/3e'0) 



In the zero size limit, a — > 0, (namely, 0o — >^ while c;, remains fixed) the 
above equation reduces to the regular Poisson-Boltzmann equation: 



VV = Cbsmh{z(5eip) (12) 



For an asymmetric l:z electrolyte the derivation is very similar and the mod- 
ified PB equation is: 



2 47rZeC5 ^zpei._^-pe^ 

^ e 1 - 00 + 00 (e^/^^^ +z e-^^^)/{z + 1) ^ 

where 0o = (2; + l)a^C6 is the combined bulk volume fraction of the positive 
and negative ions. 

Finally, if the solution is salt free and contains only negative counterions of 
valency —\z\, the modified PB equation becomes 



where 0o = a^co is the volume fraction at an arbitrary reference point ro where 
'0(ro) = and c{yq) = cq. Note that the reference point of zero potential 
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does not lie at infinity. The salt-frcc system contains only the counterions 
which neutralize the surface charges. Since the surface is taken to be infinite 
in its size, the potential does not go to zero as a; — > oo, but it diverges to 
— oo. Physically this divergence is not a problem because the electric field and 
counterion density tend to zero at large distances. 



2.2 Phenomenological Free Energy Derivation 



The modified PB equation (eqs. 11, 13 and 14) can also be derived from a 
phenomenological free energy [19]. Let us consider again the symmetric z:z 
case. This is done by expressing the free energy of the system F = Uei — TS 
in terms of the local electrostatic potential -(/^(r) and the ion concentrations 
c^{r). The electrostatic contribution is 



C/ei = / dr — — \Vi/j\'^ + zec'^ilj — zec ip — /j.+c'^ — n^c 
J I Stt 



(15) 



The first term is the self energy of the electric field and the next two terms are 
the electrostatic energies of the ions. The last two terms couple the system to 
a bulk reservoir, where ii± are the chemical potentials of the ions. 

The entropic contribution is 



T5=^ / dr 



c+a'^ ln(c+a'^) + c~a'^ ln(c"a-^) 
+ (1 - c+a^ - c-a^) ln(l - c+a^ - c-q^) 



(16) 



The first two terms represent the translational entropy of the positive and 
negative ions, whereas the last term is the entropy of the solvent molecules. It is 
this last term that is responsible for the modification of the PB equation. Note 
that it is also possible to include additional short-range (non-electrostatic) 
interaction terms in the MPB free energy [28], but we will not consider them 
here. 



Minimizing the total free energy with respect to ip and yields the Poisson 
equation 



V^-?/' = zec'^{r) — zec~{r) 



(17) 



where the ion concentrations are given by 
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^± _ '"6 ^ (18) 

1 — 00 + 00 cosh {Pzeijj) 

Here, as before, 0o = 2cba^ denotes the bulk volume fraction of the small 
ions. Combining the above two expressions recovers the modified PB equation 
(eq. 11). The same approach can also be applied in the derivation of eqs. 13, 
14 [19]. A similar expression was suggested in the 50's by Eigcn [16] and more 
recently by Kralj-Iglic and Iglic [17]. Similar ionic distributions can be also 
obtained for solid electrolytes [33,34]. 

This approach deviates significantly from the original PB equation for large 
electrostatic potentials \Peip\ 1. In particular, the ionic concentration is 
unbound in the standard PB approach, whereas here it is always bound by 
1/a^ ("close packing") as can be seen from eqs. 11, 13 and 14. This effect 
is important close to strongly charged surfaces immersed in an electrolyte 
solution. 

Note that for high positive potentials, j3eip ^ 1, the contribution of the posi- 
tive ions is negligible and the negative ion concentration follows a distribution 
reminiscent of the Fermi-Dirac distribution 



where the excluded volume interaction plays the role of the Pauli exclusion 
principle, and = is the chemical potential of the negative ions. 

The effect of the additional entropy term and the fact that it limits the concen- 
tration of ions near the surface can be demonstrated by the following simple 
argument: consider a nearly saturated region near a highly charged (positive) 
surface where c'a^ — > 1. On one hand, the system can gain entropy by pushing 
solvent molecules in between the positive ions. On the other hand, diluting 
the saturated layer costs electrostatic energy. Altogether, the excess free en- 
ergy per unit volume of diluting the saturated layer with a volume fraction 
= 1 — c^a"^ of solvent molecules is: 



A/^^r?(lnr?-l) + ^r? (20) 



il)s is the electric potential near the surface. The balance between these two 
terms results in an optimal non-zero dilution equal to 



^* ^ Q-zPe\i,s\ ^21) 

and the free energy gain due to this dilution is 
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(22) 



We conclude that entropy will always drive solvent molecules into the satu- 
rated layer, although at high surface charge densities their amount becomes 
exponentially small. 



3 Adsorption to a Flat Surface 



The modified PB equation can help to interpret the adsorption experiments of 
Cuvillier et al. [18]. In their experiments large (diameter of about lOA) nega- 
tive multivalent {z = 3 oi z = A) ions were adsorbed onto a positively charged 
Langmuir monolayer (Fig. 1). The large polyanions such as phosphotungstic 
acid (H3PW12O40) were dissolved in an aqueous subphase and attracted to a 
cationic Langmuir monolayer such as a fatty amine surfactant (C20H41-NH2), 
spread at the water /air interface. 

Consider, therefore, a solution containing an asymmetric (I'z) electrolyte in 
contact with a single planar surface of charge density a > 0. Assuming that 
the system is homogeneous in the lateral directions, the modified PB equation 
(eq. 13) reads: 



Q^y _ Q-y 

y"ix) = At:Ibzci, ■ — — (23) 

''l-(/)o + 0o(e^2/+^e-2/)/(^+l) ^ ' 

where y{x) — (3eip{x) is the reduced electrostatic potential as function of 
the distance x from the surface, and Ib — e^/eksT is the Bjerrum length. 
Ib = 7 A for an aqueous solution at room temperature. Equation 23 can be 
solved numerically yielding the electrostatic potential and ion concentration 
profiles as a function of the distance from the surface. 

Typical results are presented in Fig. 3 for various values of the ion size a 
together with the solution of the original (a = 0) PB equation. In Fig. 3a the 
negative ion concentrations are plotted and in Fig. 3b the electric potential. 
Since at high surface charge densities the positive ion concentration is small 
near the surface, only the negative ion profiles are shown. 

Clearly, the ionic concentration saturates to its maximal value in the vicinity 
of the charged surface. This should be contrasted with the original PB scheme 
which leads to extremely high and unphysical values of cj = c~(0), especially 
for multivalent ions. In the saturated region, the ionic concentration tends to 
l/a^, leading to more pronounced deviations from PB for larger ions. Note, 
that in contrast with the Stern layer model, the width of the saturated layer 
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(Fig. 3a) is not strictly equal to a. As will be demonstrated below, it depends 
not only on the parameter a but also on the surface charge density a. 

In the saturated layer the right-hand side of eq. 13 becomes a constant, and 
the electrostatic potential is quadratic 



ipyx) —ips X + 



-X 



(24) 



where ips is the surface potential and the boundary condition ip'\s — —^TTcr/e is 
satisfied. As can be seen in Fig. 3b, the parabolic curve is a good approximation 
for ■il){x) close to the surface. The width of the saturated layer I* is not strictly 
equal to a. It can be easily estimated from eq. 24 to be 



r ~ a^a/ze 
in qualitative agreement with Fig. 3a. 



(25) 



The surface potential ips can be calculated in a closed analytical form using 
the first integral of eq. 23, with a single assumption that the positive ions 
(co-ions) density is neghgible at the surface 



where 



ze 



In 



e^i -(1 - 0o)l - In(c6a3) 



ze 



Ci-ln(c6a3) (26) 



Ci 



27raV2 
eksT 



(27) 



The last approximation in eq. 26 is valid only for high values of (^i ^ 1. The 
dimensionless parameter C,i is a measure of the importance of excluded volume 
interactions. It can be written as a ratio between the volume of a single ion 
and an electrostatic volume 



a? 



Ci = ^ (28) 

where S = e/\a\ is the area per surface charge, and Aqc = S/(27r/s) ^ 
l/|cr| is the Gouy-Chapman length [35] characterizing the width of a diffusive 
electrolyte layer near a charged surface in a salt-free solution. 

Similarly, the concentration of negative ions at the surface can be calculated 
leading to a modified Grahame Equation [1] 
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(29) 



In the above equation, as in eq. 26, the only additional approximation is to 
neglect the concentration c+ of the positive (co-ions) close to the surface. At 
low surface charge -C 1, and the ion concentration reduces to the PB results 



but for high surface charge Ci ^ 1, the deviation from the PB case is substan- 
tial. 

The surface concentration of negative ions is depicted in Fig. 4a, where cj 
is plotted as a function of the surface charge density, a/e. Equation 29 is 
presented for three different ions sizes and the original PB case (eq. 30) is 
shown as well for comparison. At low surface charges the four curves are 
similar, but as the surface charge increases, the curves deviate from each other. 
The large ions deviate first and saturate to a lower surface concentration 
c~ — The deviation point corresponds to (i ~ 1. For a = lOA this gives 
cr/e ~ O.OOSA^^ in qualitative agreement with Fig. 4a. It is interesting to note 
that the ionic concentration near the surface at high surface charge density 
depends only weakly on the bulk electrolyte concentration, Cb- 

At high values of Ci, the width of the saturated layer is of the order of /* 
(eq. 25) and the amount of charge in the saturated layer can be estimated by 



zee. 



T ~ \l - {1 - (po) e-^'] a (31) 



In Fig. 4b the ratio between the ion charge density of the saturated layer a* and 
a is plotted as function of the specific surface area per unit charge, as is often 
measured in experiments. At high surface charge densities (or equivalently, 
small surface area per unit charge) the saturated layer plays a dominant role 
in neutralizing the surface charge density. As the surface charge density is 
lowered, the width of the saturated layer decreases until it vanishes. This 
occurs when is of order unity, corresponding to a* /a ~ 1 — 1/ e in the limit 
(j)Q <ti I. This crossover is indicated in Fig. 4b as an horizontal line below which 
the saturated layer is no longer well defined. It should be noted that in our 
approach the counterions never over-compensate the surface charges, namely, 
a* /a < 1. 

In the experiments of Cuvillier et al. [18] the adsorbed ion density (per unit 
area) in the solution, a*, is measured by X-ray refiectivity. It is then related 
to the surface charge density a, which is controlled by the Langmuir trough 
lateral pressure. The experiments show very clearly the presence of the steric 
effects for these large ions (of estimated size of lOA). As the surface charge 
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density increases, a* jo decreases in accord with our findings (Fig. 4b) and 
in contrast to the original PB approach. The experiments also show some 
evidence that the amount of charge in the adsorbed layer (per unit area) might 
exceed the original surface charge density and lead to over-compensation of 
surface charges [18]. Our approach does not yield such an effect. It is possible 
that additional ion- ion or ion-surface interactions are responsible for this effect. 

The two different adsorption regimes are shown in Fig. 5 where the dividing 
line Ci = 1 marks the onset of steric effects. We have seen that at low charge 
densities the parameter ^ 1 and our results coincide with those of the 
original PB equation which does not take into account the steric effects. On 
the other hand, when 3> 1, corresponding to high surface charge densities 
or large ion sizes, steric effects are, indeed, important. The concentration of 
counterions can not exceed its maximal value of close packing and this affects 
the surface potential and its dependence on the surface charge. 



4 conclusions 

In this work we have studied the effect of ion size on the density profiles 
of counterions near a charged surface. For this purpose, a modified Poisson- 
Boltzmann (MPB) equation was derived which includes the short range steric 
repulsion between ions in addition to the electrostatic interactions. A formal 
lattice gas derivation was presented as well as a phenomenological one. In the 
former derivation the steric repulsion results from the maximal occupancy of 
a single lattice site, while in the latter it stems from the translational entropy 
of the ions and solvent molecules. The two points of view are equivalent and 
result in the same MPB equation. 

The MPB formalism was applied to a simple case of relatively large coimteri- 
ons adsorbing to a strongly charged surface. This model system corresponds 
to recent experiments [18] where multivalent large ions were adsorbed to a 
charged Langmuir monolayer. The advantage of the monolayer setup is the 
possibility to control the surface charge density via application of a surface 
pressure at the water/air interface. 

We show that at high surface charge densities, in agreement with experiments, 

a saturated counterion layer is formed close to the surface. This is due to the 
fact that the ion concentration is limited by its maximal value, namely, closed 
packing. This is in contrast with the standard Poisson-Boltzmann formalism 
where the ion concentration can be arbitrarily high. 

The difference between the MPB and standard PB formalism can be also 
expressed in terms of a modified Grahame equation relating the concentration 
of counterions at the surface to the surface charge density. This comparison 
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introduces a dimensionless parameter (i = a^/SAcc where a is the ion size, 
S the specific area per surface charge and Aqc = S/27r/B the Gouy-Chapman 
length, Ib being the Bjerrum length. As long as C ^ 1 the steric efi^ects are 
weak. However, when C ^ 1 (large a, small E) they cannot be neglected any 
more and saturation appears in the adsorbed layer. 

The quadratic corrections to mean field theory presented here can be easily 
calculated. They correspond to the Random Phase Approximation (RPA). In 
a homogeneous system, it can be shown that the effective charge-charge in- 
teractions are of the Dcbye-Hiickcl type with an unmodified screening length. 
In a non-uniform system (e.g. in the presence of a charged surface) the calcu- 
lation is more involved. It requires knowledge of the analytical solution of the 
corresponding mean field equation [36]. 

Since the aim of this work was to identify and clarify the consequences of a 
single effect, namely the ionic size, we did not include additional contributions 
which might affect similar systems. Among these are specific ion-surface inter- 
actions, attractive ion- ion interactions and ion- ion correlations [22-28]. The 
latter are especially important for multivalent ions but their inclusion usually 
leads to complex integral equations where the underlying physics is not as 
transparent as in the MPB approach. 

The MPB equation can be also applied to different geometries to study the ad- 
sorption to curves surfaces [37] and the interactions between charged surfaces 
in the presence of large ions. 
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Fig. 1. Schematic view of the adsorption of large ions to a charged monolayer 
[18]. The surface charge is carried by amphiphilic molecules which arc confined to 
the air/water interface. The surface charge density can be varied continuously by 
changing the area per amphiphilic molecule. 
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Fig. 2. Schematic view of an electrolyte on a lattice model. The lattice cells are 
located at rj and assigned a spin-like variable Sj = 0, ±1. 
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Fig. 3. (a) Concentration profiles of negative multivalent ions c~ (x) near a positively 
charged surface as obtained from the numerical solution of eq. 23 for two different 
ion size a = 7.5A (dotted line) and a = lOA (dashed line). The saturated layer 
width /* ~ 2A and SA, respectively, is indicated by small arrows. The solid line 
represents the concentration profile of the standard PB equation, (b) Calculated 
electrostatic potential profiles near the surface plotted together with the parabolic 
approximation (eqs. 24, 26). The dotted, dashed and solid lines are as in (a). The 
bulk concentration is Cb = O.IM for a l:z electrolyte with z=4:. The surface charge 
density a is taken as one electron charge per 50A2. The aqueous solution with e = 80 
is at room temperature so that the Bjerrum length is Ib = 7A. 
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Fig. 4. (a) Concentration of counterions at the surface, cj = c'~(0), as a function 
of the surface charge density cr/e for different ions size, a. The PB concentration is 
also plotted for comparison, (b) Ratio of the saturated layer charge density and the 
surface charge density, a* ja = zec~l* /a^ as a function of the specific surface area 
per unit charge, e/cr, for different ion sizes. The l:z electrolyte bulk concentration is 
C5 = IrnM and the valency z = A. The horizontal line indicates the crossover region 
^1 ~ 1 in low salt concentrations. 
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Fig. 5. Schematic diagram of the adsorption regimes of the large ions to a charged 

surface. The dashed line separates between the regime where steric effects are too 
weak to be relevant and the regime where the short range ion-ion repulsion becomes 
important. The horizontal axis is the area per unit charge e/a in and the vertical 
axis is the ion size a. 
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